Immunoreactivity of LMO7 and other molecular markers as potential prognostic factors in oropharyngeal squamous cell carcinoma

Background Despite the better prognosis associated with human papillomavirus (HPV)-positive oropharyngeal squamous cell carcinoma (OPSCC), some patients experience relapse and succumb to the disease; thus, there is a need for biomarkers identifying these patients for intensified treatment. Leucine-rich repeats and immunoglobulin-like domain (LRIG) protein 1 is a negative regulator of receptor tyrosine kinase signaling and a positive prognostic factor in OPSCC. Studies indicate that LRIG1 interacts with the LIM domain 7 protein (LMO7), a stabilizer of adherence junctions. Its role in OPSCC has not been studied before. Methods A total of 145 patients diagnosed with OPSCC were enrolled. Immunohistochemical LMO7 expression and staining intensity were evaluated in the tumors and correlated with known clinical and pathological prognostic factors, such as HPV status and LRIG1, CD44, Ki67, and p53 expression. Results Our results show that high LMO7 expression is associated with significantly longer overall survival (OS) (p = 0.044). LMO7 was a positive prognostic factor for OS in univariate analysis (HR 0.515, 95% CI: 0.267–0.994, p = 0.048) but not in multivariate analysis. The LMO7 expression correlated with LRIG1 expression (p = 0.048), consistent with previous findings. Interestingly, strong LRIG1 staining intensity was an independent negative prognostic factor in the HPV-driven group of tumors (HR 2.847, 95% Cl: 1.036–7.825, p = 0.043). Conclusions We show for the first time that high LMO7 expression is a positive prognostic factor in OPSCC, and we propose that LMO7 should be further explored as a biomarker. In contrast to previous reports, LRIG1 expression was shown to be an independent negative prognostic factor in HPV-driven OPSCC.


Introduction
Oropharyngeal squamous cell carcinoma (OPSCC), most often originating from the tonsils or base of the tongue, is a relatively rare group of diseases, accounting for approximately 0.5% of all new cancer cases each year [1].However, human papillomavirus (HPV)-associated OPSCC has one of the most rapidly increasing incidence rates of all cancers [2].In the most recent TNM staging system, stratification by HPV status is recommended, in accordance with findings of distinct molecular profiles, tumor characteristics, and outcomes [2,3].The prognosis is substantially better for patients with HPV-driven tumors, even though approximately 15% of them have recurrent disease [4].HPV status does not affect current treatment recommendations.Still, considering the younger age and better prognosis associated with HPV-driven tumors, means are taken to de-escalate treatment to increase the quality of life for this group of patients [5].Identifying HPV-positive patients with a more aggressive disease remains a challenge, and further prognostic and predictive markers are thus needed.
The human LRIG (leucine-rich repeats and immunoglobulin-like domains) family of integral transmembrane proteins comprises LRIG1, LRIG2, and LRIG3 [6].LRIG proteins are ubiquitously expressed in human tissues and even though their role in tumorigenesis is not fully understood, LRIG proteins have a prognostic impact in various tumor diseases (reviewed in [7,8]).LRIG1 functions as a tumor suppressant and has proven to be a negative regulator of receptor tyrosine kinase signaling [9], as well as a promotor of bone morphogenic protein (BMP) signaling [10].LRIG1 expression has been shown to correlate with HPV status in OPSCC [11] as well as in other cancers [12,13].In HPV-driven OPSCC, LRIG1 expression was also proven to be an independent prognostic marker associated with longer survival [11].
Studies have shown that there is an interaction between LRIG1 and the LIM domain 7 protein (LMO7) [14].LMO7 has a stabilizing role in adherence junctions as an afadin-and a-actinin-binding protein [15], is involved in cell migration [16], and acts as a transcription factor of muscle-related genes [17].It has limited expression in normal tissues but is expressed in different cancers and is suggested to have an important function in cancer development and metastasis [16,18,19].Its expression as a prognostic marker has been studied in several tumor diseases, with differing results.In LRIG1-positive non-small cell lung cancer, LMO7 was a negative prognostic marker [14], while another study of lung adenocarcinoma showed that the loss of LMO7 was associated with a worse prognosis [18].It was also deemed tumorpromoting in pancreatic cancer [20].In vulvar cancer, high LMO7 expression was associated with better overall survival (OS) and disease-free survival (DFS) [21].This connection was stronger in non-HPV-driven tumors [21].LMO7 has not been studied in OPSCC, and thus, the aim of this study was to elucidate the role of LMO7 in OPSCC and investigate its association with LRIG1 and other known prognostic clinicopathologic characteristics.

Patients and specimens
Patients undergoing treatment for oropharyngeal cancer at Umeå University Hospital between 1990 and 2013 were identified via the Swedish Cancer Registry and patient records.Formalin-fixed, paraffin-embedded (FFPE) tumor tissue samples were retrieved from the biobank at Västerbotten County Council after approval from the Regional Ethical Review Board at Umeå University (Dnr 2013-475-31 M).Informed consent from participants alive at the start of the study was obtained.Clinical data were collected from patients' records by a medical doctor.Out of initially 329 eligible patients treated with an intention to cure, 36 declined to participate, 34 were excluded due to insufficient clinical data and in 114 cases material could not be obtained.Thus, 145 patients were included in the study.Individual patients may be excluded from the respective molecular analyses due to missing data/ lack of material.

HPV DNA analysis
DNA was extracted from FFPE sections using DNA Tissue FFPE Tissue Kit or QIAamp MINI Kit (Qiagen, Inc., Valencia, CA, USA) according to the manufacturer's instructions.DNA samples were stored at -20°C.HPV-PCR analysis was carried out using 100 ng of extracted DNA from each patient.General primers GP5 + /6 + were used to amplify a fragment of the conserved HPV L1 gene as previously described [22].Amplification was performed in a Biometra Professional Thermocycler (Thermo Fisher Scientific, Waltham, MA, USA) or a T100 Thermal Cycler (Bio-Rad, Hercules, CA, USA).PCR products were run on a 2.5% agarose gel in a buffer consisting of 50 mM Tris/37 mM Borate/1.3mM EDTA and stained with 0.5X GelRed (Biotium, Hayward, CA, USA).Gels were visualized under UV light.Fragments of 130-150 bp were considered HPV-positive.To avoid false HPV-negative results due to disruption of the L1 gene when using the GP5 + /6 + primers, negative samples were reanalyzed with the general primers CpI/IIG as described [23].Samples with products of approximately 188 bp were then considered positive.
HPV typing using HPV16 type-specific primers was also performed, and samples positive for HPV with general primers and negative for HPV16 were sequenced for determination of HPV type on the amplicon generated from these primers, as previously described [24].

Evaluation and classification of immunostaining
LRIG1, LMO7, CD44, and p53 staining was scored based on both staining intensity and percentage of immunoreactive epithelial cells.The intensity was evaluated on a four-grade semiquantitative scale; 'no staining' , 'weak' , 'intermediate' or 'strong intensity' , and the fraction of positive cells was scored separately and evaluated as 0%, 1-25%, 26-50%, 51-75%, and 76-100%.Ki67 staining was scored as the percentage of tumor cells with positive nuclear staining and the tumors were subgrouped into 'at or above mean' and 'below mean' .The fraction and intensity of p16-positive cells were evaluated and only cases with strong and continuous nuclear and cytoplasmic expression in > 70% of positive cells were considered p16-positive [25].Cases with cytoplasmic staining only were considered negative.Tumors with no or high (> 75%) P53 expression, or strong staining intensity, were considered p53 mutated in accordance with a previous study [26].All evaluations were performed by an independent pathologist, with experience in interpreting LMO7 and LRIG staining, blinded to the clinical data and outcomes.

Statistical analyses
DFS was defined from the date of completed treatment to the date of recurrence (according to pathological anatomical diagnosis or date of first recorded clinical progression).Death without documented recurrence was censored at the date of death.OS was defined as the time from the date of completed treatment to the date of death irrespective of the cause.If a patient was still alive at the time of accessing the patient's records, the case was censored at that date.
Patient characteristics were analysed with the independent samples t-test or one-way ANOVA for comparison of means, or the Pearson chi-square test or Fisher's exact test for ordinal variables.Two-sided P values were reported.DFS and OS were illustrated in Kaplan-Meier graphs and a log-rank test was used for comparison of the survival probabilities.Univariate Cox proportional hazard models were used to calculate the hazard ratio (HR).All variables with significant results in the univariate analysis were further tested in a Cox regression multivariable analysis to evaluate the independent factors' influence on the risk of death.To avoid overfitting, in some cases, variables with significant results in the univariate analysis were further investigated through backward stepwise elimination to estimate adjusted HRs.P values < 0.05 were considered statistically significant.All statistical analyses were performed using SPSS 28 software (IBM SPSS Statistics 28, IBM, Armonk, NY, USA).

Study population
A total of 145 patients treated for oropharyngeal cancer with curative intent at Umeå University Hospital, Sweden, between 1990 and 2013 from whom clinical data and representative material could be obtained were enrolled in the study.The main characteristics of the study patients are summarized in Table 1.The patients were initially staged according to the 7th edition of the TNM staging classification and later restaged according to the Union for International Cancer Control (UICC) 8th edition [3].p16 positivity causes downstaging of the tumor compared to the previous staging system, and this was also the case in our study (Table 1).The nodal status of the patients was as follows: N0 = 28, N1 = 11, N2 = 81, N2a = 2, N2b = 5, N2c = 6, N3 = 9, and N3b = 1.All patients were treated with radiotherapy (conventional or accelerated); in two cases the treatment was interrupted, and the intention was changed to palliation.Six patients underwent diagnostic surgery before radiotherapy, and 54 patients underwent surgery (neck dissection), the majority following radiotherapy.One patient had neoadjuvant chemotherapy, and nine patients were treated with concomitant chemotherapy.A total of 14/143 (10%) patients had residual tumor at the first follow-up, and 34/143 (24%) relapsed during the study period.There were 61 deaths from any cause, of which 47 patients died from OPSCC (Table 1).

HPV and p16 analyses and patient characteristics
Seventy-seven percent (111/143) of the tumors were HPV-positive, 81% (117/143) were p16-positive and 76% (108/143) were positive for both HPV and p16.Two tumors were HPV-positive and p16-negative.Eight tumors were p16-positive but HPV-negative.In two cases, HPV data could not be retrieved, and in two cases, data on p16 were missing.In two cases one variable was positive, but the other was missing; thus, these cases could not be categorized as HPV-driven or not and excluded from further analyses.p16 is used as a surrogate marker for HPV infection and in the UICC-8 staging system for the classification of OPSCC [3].The p16 analysis was thus used for restaging.However, in this study, tumors positive for both HPV and p16 were regarded as HPV-driven, in accordance with [27].From the 111 HPVpositive tumors, genotyping was performed on 92.In 19 cases, sufficient material was lacking.HPV16 was the most frequent type, found in 84 cases (75%).HPV33 was found in four cases, and HPV18, HPV26, HPV35, and HPV45 were found in one case each.HPV-driven tumors were significantly associated with earlier stages, tonsillar origin, lower differentiation grade, and lower percentage of smokers compared to the patients with non-HPVdriven tumors, and all HPV-driven tumors except one exhibited p53 wild-type staining pattern (Table 1).
There was no significant difference in the expression between HPV-and non-HPV-driven tumors (Table 1).
When the staining was dichotomized as no/weak versus intermediate/strong, no/weak/intermediate versus strong, and 0-50% versus 51-100% no significant difference was found.LMO7 expression correlated with LRIG1 expression (p = 0.048), differentiation grade (p = 0.036) and vital status (p = 0.022) (Table 2).LMO7 expression did not correlate with any other molecular or clinical feature (Table 2).Survival analyses showed that a high (51-100%) staining percentage of LMO7 was associated with favorable OS (p = 0.044) (Fig. 1c).Stratification into HPV-/non-HPV-driven tumors did not reveal any difference within any of the groups.To further explore these results, we divided the patients into four groups according to their combined HPV and LMO7 status; HPVdriven tumors with high (51-100%) LMO7 expression (HPV + , LMO7 +), non-HPV-driven tumors with low (0-50%) LMO7 expression (HPV-, LMO7-), HPV-driven tumors with low LMO7 expression (HPV + , LMO7-) and non-HPV-driven tumors with high LMO7 expression (HPV-, LMO7 +).The clinical characteristics of the different groups can be viewed in Table 3. Survival analysis revealed that HPV + , LMO7 + tumors were associated with the best, and HPV-, LMO7-tumors with the worst, OS (Fig. 3a, p = < 0.001).Patients with HPV-driven tumors that expressed a low percentage of LMO7 had the second-best survival and non-HPV-driven tumors with high LMO7 expression the third best.Analysis of DFS gave the same results (Fig. 3b, p = < 0.001).There was no difference in OS or DFS between the tumors when grouped according to no/weak versus intermediate/ strong staining intensity.In the cohort of 143 patients, immunohistochemical LRIG1 expression was evaluated by staining intensity (exemplified in Fig. 2d-f ) in 96 tumors and the fraction of positive cells in 94 tumors.LRIG1 staining was negative in 1 (1%) tumor, and the staining intensity was weak in 4 (4%), intermediate in 75 (78%), and strong in 16 (17%); 2 (2%) had 26-50% positive cells, 5 (5%) had 51-75% positive cells, and 86 (92%) had 76-100% positive cells.There was no significant difference in the expression between HPV-and non-HPV-driven tumors (Table 1).When the staining was dichotomized as absent/weak/intermediate versus strong or the percentage of positive cells as 0-75% versus 76-100% there was no significant difference in the expression between HPV-and non-HPV-driven tumors.
There was no difference in OS or DFS when comparing tumors with a strong LRIG1 staining intensity to no/ weak/intermediate.However, using HPV status as strata revealed a significantly worse OS in the patients with HPV-driven tumors when the tumors presented with a strong LRIG1 staining intensity (p = 0.016) (Fig. 1d).The patients were further grouped according to their combined HPV and LRIG1 status; HPV-driven tumors with strong LRIG1 staining intensity (HPV + , LRIG1 +), non-HPV-driven tumors with no/weak/intermediate LRIG1 staining intensity (HPV-, LRIG1-), HPV-driven tumors with no/weak/intermediate LRIG1 staining intensity (HPV + , LRIG1-) and non-HPV-driven tumors with strong LRIG1 staining intensity (HPV-, LRIG +).The  4. Survival analysis showed that the HPV-driven, LRIG1 staining intensity no/weak/intermediate tumors were associated with the best OS (Fig. 3c, p = < 0.001).The group with HPV-driven tumors with strong LRIG1 staining intensity had a better OS than the two groups of non-HPV-driven tumors.Studying DFS, the HPV-driven tumors had a beneficial prognosis, regardless of LRIG1 status (Fig. 3d, p = 0.003).When the percentage of positive cells was dichotomized as 0-75% versus 76-100% there was no significant difference in OS or DFS.No further difference was seen within the group of HPV-or non-HPV-driven tumors.
Ninety-five tumor tissue samples were successfully stained for CD44 and evaluated both by staining intensity and fraction of positive cells (as exemplified in Fig. 4ac) (Table 1).When the intensity was dichotomized as absent/weak versus medium/strong, the HPV-driven tumors more often expressed a lower staining intensity (33/70, 47% vs. 6/26, 24%, p = 0.043).When the fraction of positive cells was dichotomized as 0-25% versus > 25%, the HPV-driven tumors also expressed CD44 to a significantly lower extent (35/70, 50% vs. 21/25, 84%, p = 0.003).There was a significant correlation between strong CD44 staining intensity and poor overall survival (Fig. 5a, p = 0.016), which was also observed when analyzing the HPV-driven group of tumors separately (Fig. 5b, p = 0.02).Patients with tumors exhibiting a strong CD44 staining intensity also had a significantly shorter DFS (Fig. 5c, p = 0.011).There was no difference when CD44 intensity was dichotomized as no/weak versus intermediate/  strong or the percentage of positive cells as 0-25% versus 26-100% or 0-50% versus 51-100%.Strata according to HPV status did not reveal further significant results regarding OS or DFS.The expression of Ki67 was evaluated in 96/143 tumors (exemplified in Fig. 4d-f ).The mean fraction of cells expressing Ki67 was 52.5%.In the HPV-driven tumors the mean was 56% (n = 70) and in the non-HPV-driven tumors the mean was 43% (n = 26).The HPV-driven tumors expressed Ki67 to a significantly higher extent than the non-HPV-driven tumors (Table 1).Using the mean as a cutoff, a high expression of Ki67-positive cells correlated with longer OS (Fig. 5d, p = 0.024).Stratifying patients according to HPV status gave the same results in the patient group of non-HPV-driven tumors (Fig. 5e, p = 0.04).No significant results were observed when analyzing DFS.
P53 expression intensity and the number of positive cells were evaluated in 97/143 tumors (exemplified in Fig. 4g-i).There was a significant difference between HPV-driven and non-HPV-driven tumors.All HPVdriven tumors except one (1.5%),exhibited p53 wildtype pattern, in contrast to the non-HPV-driven tumors where 18 (69%) were p53 mutated (Table 1).Survival analyses showed that p53 wild-type expression correlated with significantly longer OS (Fig. 5f, p = < 0.001) and DFS Vital status Alive 30 (52.5) 1 (Fig. 5g, p = < 0.001), in contrast to mutated p53 expression.Using HPV status as strata revealed that p53 wildtype expression was associated with longer OS (Fig. 5h, p = 0.021) also in the subgroup of non-HPV-driven tumors.No survival analysis was performed in the group of HPV-driven tumors since there was only one case of p53 mutation.

Univariate analyses
The crude hazard ratio for overall survival is shown in Table 5.There was no association between OS and sex, smoking, or LRIG1 staining intensity.There was a significant association between age, tumor localization, differentiation grade, stage according to TNM 8, HPV status, LMO7 expression, CD44 staining intensity, Ki67 and p53 expression, and OS.Patients with a high percentage (51-100%) of LMO7-positive cells had a favorable OS (HR 0.515, 95% CI: 0.267-0.994,p = 0.048).The univariate associations of patient and disease characteristics with disease-free survival are shown in Table 6.Stage according to TNM 8, HPV status, CD44 staining intensity, and p53 expression correlated with DFS.Grouping the patients according to HPV status revealed that in the HPV-driven tumors (108/143) (Table 7), there was an association between OS and age, stage according to TNM 8, differentiation grade, LRIG1, and CD44 staining intensity.There were no significant associations between DFS and the different patient and disease characteristics (data not shown).In the non-HPV-driven tumors (35/143), univariate analyses showed that p53 mutated tumors were associated with shorter OS (HR 3.881, 95% Cl: 1.124-13.404,p = 0.032, data not shown).There were no other significant association with DFS or OS (data not shown).

Multivariate analysis
Parameters with a significant correlation to OS or DFS, respectively, in the univariate analyses were included in the multivariate analysis.To avoid overfitting regarding overall survival, backward stepwise selection was used to select which of the variables to include in the final multivariate analysis (Tables 5 and 7).This analysis revealed that stage according to the 8th edition of the UICC TNM staging system and Ki67 expression were the only independent prognostic indicators of OS.
There were no independent prognostic factors for DFS (Table 6).In the group of patients with HPV-driven tumors (Table 7), the multivariate analysis revealed that age, stage according to TNM 8, and LRIG1 staining intensity were independent prognostic factors of OS.Strong LRIG1 staining intensity was associated with a worse prognosis (HR 2.847, 95% Cl: 1.036-7.825,p = 0.043).four groups with significantly different OS (Fig. 3a).The beneficial effect of a high LMO7 expression may in part explain the difference in prognosis seen within the group of patients with HPV-driven tumor.When grouping was performed by combining HPV and LRIG1 status the difference was not as pronounced.We thus suggest that LMO7 primarily needs to be further explored as a potential biomarker for risk-stratifying patients within the HPV-driven and non-HPV-driven OPSCC tumor groups respectively.LMO7 was a positive prognostic factor for OS in OPSCC in univariate analysis (HR 0.515, 95% CI: 0.267-0.994,p = 0.048).LMO7 expression also correlated with LRIG1 expression (p = 0.048), consistent with previous findings [14].Interestingly, LRIG1 was an independent negative prognostic factor in the HPV-driven group of tumors (HR 2.847, 95% Cl: 1.036-7.825,p = 0.043).
No conclusions on LMO7 expression could be drawn when stratifying for HPV status.In the group of non-HPV-driven tumors, p53 expression was the only significant prognostic factor.However, because of the limited number of cases (n = 35) in the analyzed cohort, the possible role of LMO7 as a prognostic marker in this group cannot be ruled out.The connection between LMO7 and HPV has to our knowledge previously only been studied in vulvar cancer [21].The authors found that a high LMO7 expression was associated with a favorable OS in the group of HPV-negative patients (n = 70).Considering this and our results from the present study, the connection between HPV and LMO7 should be further investigated in a larger cohort of OPSCC patients.
Table 7 Cox regression univariate and multivariate associations of overall survival with clinicopathologic characteristics, in HPV-driven tumors LRIG1 is generally considered a positive prognostic factor, as seen in vaginal carcinoma [12], prostate [28], breast [29,30], lung [31], and uterine cervical cancer [13,32], among others [7].LRIG1 is also a positive prognostic marker in OPSCC; Lindquist et al. (2014) found that LRIG1 was associated with longer OS and DFS in OPSCC in the whole cohort and in the subgroup of HPVdriven tumors [11].In contrast, our results indicate that LRIG1 staining intensity is a negative prognostic factor in HPV-driven tumors.The LRIG1 antibody used in this study has previously been used in an investigation of LRIG1 expression in vulvar cancer.In line with our results, this study showed that no/weak LRIG1 staining in HPV-negative tumors was associated with a better OS [21].The discrepancy between these two reports and previous findings might be due to methodological reasons.Both the LRIG1-Vina antibodies used in this cohort and the LRIG1-151 antibodies used in the study of Lindquist et al. were polyclonal, targeting the same antigen, and the same procedure was used in both studies.Thus, we can justify the result as follows: evaluation of staining intensity is highly subjective, and some factors influence the staining, such as the duration of storage of the tumor section, fixation method of the tumor tissue, and methodology of the immunohistochemistry.Our results indicate that there is a need for more standardized antibody development and interpretation of the result, which needs to be tested in a larger material before conclusions can be drawn on LRIG1 in OPSCC to be used in a clinical setting.Nonetheless, high LRIG1 expression is not consistently favorable and it has been suggested that its function can be either tumor-promoting or suppressive depending on the cellular context [33].In an untreated prostate cancer cohort, LRIG1 was a marker for poor survival [28].In contrast, LRIG1 was a positive prognostic factor following prostatectomy, further indicating that the association between survival and LRIG1 is contextdependent [28].Blood levels of LRIG1 have only recently been studied, and interestingly, the results from a highgrade serous ovarian cancer cohort showed that high levels were associated with poor prognosis [34].The authors speculate that this could be connected to the BMP-promoting effect of LRIG1, as BMPs are associated with tumor-stimulating functions in ovarian carcinoma [34].Neither the role of BMPs nor the association between BMPs and LRIG1 has been studied in OPSCC.Overall, the conflicting results on the role of LRIG1 in OPSCC warrant further studies.
As expected, HPV-driven tumors were associated with longer DFS and OS (p = < 0.001).The current staging system stratifies patients according to HPV status based on p16 expression via IHC, which is used as a surrogate marker [35].However, approximately 20% of p16-positive tumors may lack transcriptionally active HPV, which may affect outcome [36].It has also been proposed that transient HPV infection may cause cases where HPV DNA can be detected but p16 is negative [37].Studies indicate that a positive assessment of p16 (via IHC) should be confirmed by high-risk HPV DNA detection [36], and this is recommended in the current Swedish guidelines [38].Thus, we chose to consider p16-positive/HPV-positive tumors as HPV-driven.
CD44 is a transmembrane glycoprotein overexpressed in several tumors.Binding of its main ligand hyaluronan causes proliferation, survival and motility-enhancing signaling [39].Absent/weak CD44 staining intensity has been found to correlate with a longer OS and DFS in OPSCC [40,41], and this was also confirmed in our material (Fig. 5a-c), as was the finding that HPV-driven tumors express both a lower CD44 staining intensity and a lower percentage of positive cells than non-HPV-driven tumors (Table 1) [40].The HPV-driven tumors expressed Ki67 to a significantly higher extent than the non-HPVdriven tumors (Table 1), consistent with previous studies [42].Surprisingly, our data suggest that Ki67 expression at or above the mean is associated with favorable prognosis, also in multivariate analysis (Table 5), in contrast to other studies [42,43].There is no consensus on what cutoff to use which makes it difficult to compare the results and thus this needs to be further investigated.Different studies use different cutoffs regarding p53 staining intensity and expression [44,45].Two reports on HPV-associated gynecological cancers suggest that tumors with no p53 staining, strong staining intensity, or a high expression of positive cells should be considered p53 mutated [46,47].This division has also been used in a recent report on oropharyngeal carcinoma [26].In concordance with this study [26], all HPV-driven tumors, except one, presented with p53 wild-type expression (Table 1).Survival curves illustrated worse OS and DFS in tumors with mutated p53 expression (Fig. 5f-g), and univariate Cox regression revealed a significant association between mutated p53 expression and worse OS (Table 5) and DFS (Table 6), as seen before [44].

Conclusions
LMO7 immunoreactivity was found for the first time to be a positive predictor of OS in OPSCC, and its role as a biomarker should be further investigated.In the group of HPV-driven tumors, strong LRIG1 staining intensity was associated with worse OS and it was shown to be a negative prognostic factor in univariate and multivariate analysis, in contrast to a previous report [11].Further studies are needed on the mechanisms and clinical applicability of LMO7 and LRIG proteins in OPSCC.

Table 1
Correlations between HPV status and different clinicopathologic characteristics

Table 1
(continued) HPV human papillomavirus, NK not known, LRIG leucine-rich repeats and immunoglobulin-like domain, LMO7 LIM domain 7 protein a Independent samples t-test b Chi 2 -test c Fisher's exact test

Table 2
Correlations between LMO7 expression and staining intensity and clinicopathologic characteristics

Table 3
Correlations between HPV and LMO7 status and clinicopathologic characteristics

Table 6
Cox regression univariate and multivariate associations of clinicopathologic characteristics with disease-free survivalHPV human papillomavirus, LMO7 LIM domain 7 protein, LRIG leucine-rich repeats and immunoglobulin-like domain * = P < 0.05